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SUMMARY 


The  manufacturing  of  co-extruded  grains  of  highly  filled  propellant  suspensions  is  a 
complicated  operation,  which  requires  a  detailed  and  realistic  understanding  of  the  various  types 
of  interface,  surface,  and  bulk  instabilities.  This  understanding  is  necessary  to  be  able  to 
successfully  manufacture  fast  core  propellants  using  co-extrusion  (ram  or  twin  screw  extrusion 
based).  In  this  investigation,  the  fundamentals  of  the  co-extrusion  process  were  investigated  by 
a  combination  of  experimental  and  simulation  studies.  The  project  included  the  development  of 
a  novel  experimental  apparatus,  which  was  used  to  investigate  various  types  of  flow  instabilities 
and  to  collect  basic  data  on  various  aspects  of  the  fundamentals  of  the  co-extrusion  process. 
Furthermore,  an  analytical  model  of  the  process  was  developed  by  Prof.  Dilhan  M.  Kalyon  of  the 
Stevens  Institute  of  Technology  (SIT),  Hoboken,  New  Jersey  to  allow  the  determination  of  the 
pressure  drop  versus  flow  rate  relationships  along  with  the  velocity  and  shear  stress  distribu¬ 
tions  when  two  fluids  are  co-extruded  side  by  side.  The  viscoplasticity  and  the  wall  slip  of  the 
two  suspensions  could  be  accommodated.  The  agreement  between  the  theory  and  the 
experimental  results  was  acceptable.  In  this  final  report,  the  techniques  and  the  results  obtained 
are  outlined  along  with  a  listing  of  the  reports  and  presentations  (refs.  1  through  8)  that  were 
made  and  a  list  of  the  lessons  learned. 


BACKGROUND  AND  OBJECTIVES 

The  manufacturing  of  co-extruded  grains  of  highly  filled  propellants  is  a  complicated 
operation,  which  requires  a  realistic  understanding  of  the  various  types  of  interface,  surface,  and 
bulk  instabilities.  This  understanding  is  necessary  to  be  able  to  successfully  manufacture  fast 
core  propellants  using  co-extrusion  (ram  ortwin  screw  extrusion  based).  In  this  investigation,  a 
novel  experimental  apparatus  was  developed  and  used  to  investigate  various  types  of  flow 
instabilities  and  to  collect  basic  data  on  various  aspects  of  the  fundamentals  of  the  co-extrusion 
process.  Furthermore,  an  analytical  model  of  the  process  was  developed  by  Prof.  Kalyon  to 
allow  the  determination  of  the  pressure  drop  versus  flow  rate  relationships  along  with  the 
velocity  and  shear  stress  distributions  when  two  fluids  are  co-extruded  side  by  side.  The 
viscoplasticity  and  the  wall  slip  of  the  two  suspensions  could  be  accommodated.  In  this  final 
report,  the  techniques  and  the  results  obtained  are  outlined  along  with  a  listing  of  the  reports 
and  presentations  that  were  made. 


EXPERIMENTAL  SET-UP 

A  modified  ram  extrusion  process  was  used  to  generate  various  types  of  co-extruded 
grains  using  polydimethyl  siloxane  (PDMS)  and  suspensions  of  PDMS  with  hollow  glass 
spheres,  with  the  experiments  designed  to  probe  the  different  types  of  instabilities,  which 
develop  during  the  co-extrusion  process  (figs.  1  to  6).  Various  typical  results  obtained  with 
differing  combinations  of  fluids  (the  combination  of  pure  binder  on  one  side  and  a  suspension  on 
the  other  side  is  the  most  difficult  and  combinations  of  suspensions  are  easier  to  handle  since 
the  elasticity  and  shear  viscosity  differences  are  not  as  pronounced  as  in  pure  binder  on  one 
side  of  the  co-extrusion  die)  are  shown  in  figures  7  to  18.  These  results  were  elucidated  in  detail 
in  our  presentations  and  progress  reports  (refs.  1  through  8)  and  thus  will  not  be  detailed  here. 
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Conclusions  from  the  Experimental  Studies 

The  experimental  results  outlined  in  our  earlier  results  indicate  that  differences  in  the 
shear  viscosity  of  the  different  phases,  which  are  used  in  the  co-extrusion  process,  affect  the 
development  of  various  types  of  interface  and  surface  irregularities.  If  the  shear  viscosity 
material  functions  of  the  two  phases  are  close  to  each  other,  then  the  effects  of  various  types  of 
instabilities  are  minimized. 

On  the  other  hand,  if  the  differences  in  the  shear  viscosity  material  functions  of  the  two 
phases  are  significant,  then  the  ram  extrusion  process  generates  a  time  dependent  change  in 
the  interface  configuration  along  with  differences  in  the  surface  irregularities  of  the  two  sides  of 
the  rectangular  extrudates.  Thus,  one  of  the  major  challenges  of  the  selection  of  the  different 
phases  to  be  used  in  the  fast  core  propellants  will  be  to  match  the  shear  viscosity  material 
functions  of  the  fast  and  slow  burn  propellants  by  correct  selection  of  temperatures  and  shear 
rate  distributions  in  the  co-extrusion  hardware.  Some  of  the  important  issues  in  the  design  of  the 
hardware  are  also  identified. 


MATHEMATICAL  MODEL  OF  THE  CO-EXTRUSION  PROCESS 

An  analytical  mathematical  model  of  the  co-extrusion  of  two  suspensions  flowing  side  by 
side  in  a  co-extrusion  die,  subject  to  the  viscoplasticity  of  the  suspensions  and  their  slip  at  the 
walls  of  the  co-extrusion  die  was  derived  by  Prof.  Kalyon.  The  basic  tenets  of  this  model  are 
described  next: 

•  There  are  two  suspensions  that  are  viscoplastic  which  are  flowing  side  by  side 
in  a  slit  die  (fig.  19).  This  is  the  same  geometry  that  is  used  in  our  experiments 
and  is  the  most  common  general  configuration,  the  modifications  of  which  will 
provide  the  other  configurations  including  the  slow/fast/slow  configuration. 

•  Both  of  the  suspensions  are  subject  to  wall  slip  (fig.  20)  at  the  walls  of  the  co¬ 
extrusion  die. 

•  The  flow  is  isothermal  and  fully-developed  (fig.  21)  and  contains  multiple 
zones  of  deforming  and  non-deforming  zones  (figs.  22  through  24).  Since  the 
constitutive  equation  differs  on  shear  stress  less  than  the  yield  stress  and  the 
shear  stress  greater  than  the  yield  stress,  and  since  the  case  of  increasing 
velocity  as  the  distance  increases  versus  the  case  of  decreasing  velocity  as 
the  distance  increases  need  to  be  treated  separately  multiple  cases  need  to 
be  addressed  (fig.  25). 

•  The  two  fluids  are  incompressible. 

The  solutions  indicate  that  there  are  twelve  cases  that  are  possible  and  each  of  these 
cases  requires  a  separate  solution.  The  simplifications  of  the  model  presented  here  would  also 
include  as  a  subset  the  flow  of  a  single  suspension  subject  to  viscoplasticity  and  wall  slip  in  a  slit 
die.  The  12  cases  are  shown  in  figure  25. 
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In  the  following  the  equations  of  the  velocity  distributions,  shear  stress  distribution,  the 
nomenclature  and  the  differentiators  between  the  twelve  cases  are  provided  for  future 
documentation  and  the  use  of  ARDEC.  The  derivation  is  novel  and  will  be  published  with 
ARDEC  acknowledgement  and  permission  in  the  future. 


Nomenclature 


*bi 

*02 

mj 

m2 

«1 

n2 

s2 

A 

Pi 

yt 


Subscript 

y 


Superscript 

I 

II 

III 

IV 

V 


=  Yield  stress  value  of  fluid  1 
=  Yield  stress  value  of  fluid  2 
=  the  consistency  index,  fluid  1 
=  the  consistency  index,  fluid  2 
=  power  law  index,  fluid  1 
=  power  law  index  of  fluid  2 
=  1/rtj  is  the  reciprocal  power-law  index  of  fluid  1 
=  \/n2  is  the  reciprocal  power-law  index  of  fluid  2 
=  Navier's  slip  coefficient  for  fluid  1 
=  Navier's  slip  coefficient  for  fluid2 
=  distance  of  the  interface  from  the  wall 


=  y-direction 

=  z-direction  (flow  direction) 


=  deformation  zone  for  fluid  1 
=  plug  flow  zone  for  fluid  1 
=  first  deformation  zone  for  fluid  2 
=  plug  flow  zone  for  fluid2 
=  second  deformation  zone  for  fluid  2 

Equations  used: 


*yz  =±r0-m 


dV \ 

"~Y  dVz] 

dy 

{  dy  J 

A,  =- 


dp  1 
dz  m\ 


5]  + 1 


Ah  —  —I 


dp  1 

dz  m2  J 


^2  + 1 
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The  assumptions  made  for  these  derivations  are: 

•  Fully  developed  flow 

•  Isothermal  conditions 

•  The  wall  slip  at  both  the  wall  is  different 

•  The  fluid  is  a  viscoplastic  and  is  represented  by  the  Herschel-Bulkley  three 
parameter  model 

•  The  fluid  is  incompressible  fluid 

•  There  is  only  laminar  flow 

•  The  more  viscous  fluid  is  at  the  bottom  called  fluid  1  and  the  less  viscous  fluid 
is  at  he  top  and  is  called  fluid  2 

A  Herschel-Bulkley  type  viscoplasticity  with  three  parameters,  i.e.,  the  yield  stress,  r0 ,  the 


consistency  index,  m  and  the  shear  rate  sensitivity  index,  n ,  i.e.,  xyz  =  ±%  - m. 


dVl, 


dy 


11:  - 1 


dV‘z 


for  TyZ 


{  dy  . 
(dVlz/dy)= 0  for 


>  x0  (-  sign  is  used  for  negative  shear  stress,  xyz ,  and  the  shear  rate 
|r>>2|  <  tq  where  i  represents  different  zones  in  the  fluids. 


There  are  12  cases.  In  the  following  the  details  of  the  derivation  of  the  first  case  is 
provided  along  with  the  final  results  for  the  other  1 1  cases.  The  distinguishing  features  of  the  12 
cases  are  also  included. 

Case  I  -  Details  of  the  Derivation  (All  Zones) 

Condition 


o<4  <yi 


yi  2<^i-h 

Zone  1 


d*y^_  _  dp 
dy  dz 


on  integration,  we  get 


I  dp  i 

=  y  +  c\ 

dz 


yz 
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iv  dp 


■*  =-ny+Cl 


IV 


v  dp  V 

**  =~T,y+c' 


From  the  continuity  of  shear  stress  c{  =  c(  =  c{  =  c'v  =  c\ ,  specifying  ryz  at  one  location  is 
adequate  at  y  =  A  ,  xyz  =  0 


dz 


Case  1  -  Diagram  (all  five  zones) 
In  the  previous  case 


(0)1  >  To,: 

i  at 

y-yb 

Tyzd) 

At  y  = 

>  l7yz| = 

:  _T0,I 

aty  = 

^2- 

M= 

-t0,2 

aty  = 

'h  > 

W= 

t0,2 

at  y  = 

0,  | 

T«|  =  T, 

,(0) 

aty  = 

h,\ 

M  =  - 

-TyZ(h) 

r>’z(  2) 


at  x=y=h±k 
2 


yz 


=  0 


Note:  0</^<y(-  and  yf  </t 


5 


We  know  that  ryz 


aty  =  Ai,ryz  =-t0(i. 


Hence,  ^  = 


-*o,i  ~^(Q) 
(-d/?/dz) 


Note:  —  <0 
dz 

Similarly, 

^3  = 

Va)=“X/,+t>-z(0) 

We  need  to  determine  one  of  these  parameters.  This  can  be  done  by  finding  the  velocity 
profiles. 

1.  VlPlug=M-TyZ( 0ifb>  +A,^I+I  at  y  =  0 

2-  {yi)--^2(^2~ yiY1+]  +vz,piug  at  y.  < y < Aj 

^z!plug  =  A  (“  Tyz  (O))^'  +  A,/^1  +A2(/^2  -yi)''2+1 

At  y  =  /t  V"=A2(^(# 

vz!plug  =  Pi iTyz(hjfb2  +  A2(h- 

and  at  y  =: 

Vz!plug  =  +Ai^l+I  +  A2(/i2  -J',')'2  +  l 

Tyz(0)  is  determined  by  equating  the  two  previous  equations  to  give 

/hfayzihifh*  +A2(/i-^3)'?2+1  =  fi(-Tyi(0)Yb'  +A|/^I+1  +  A2(/^  -  _y,)J2+l 
From  known  t>iz(0)  ,  we  determined  A\ ,  ,  and  X$  from  the  previous  equations. 


(-dpi  dz) 

-*0,2  -y(0) 

(-dp!  dz) 
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Zone  I 


0<  y <  A\ 


v/=A,L^+1-(^-^,+1J-/?iV(0) 


Zone  II 

A^^y^yi 

v/^Aj^-AVO) 


Zone  III 

yi<y<A2 


VZH  =A1^'+I -P\t, 


(0)  +  A2L(/i2 -y()'y2+1  -(A?  -  y) 


yz 


Zone  IV 

>^2  <  V  <  >^3 

VZV  =A1/i[l+l  -A(Tyz(0))+A2(y^2  _);t)'52+l 


Zone  V 

Aj<  y<h 

V/  =  A,^1  +'  -  Ar,,  (0)  +  A2L(^2  -  -  (j  -  % +‘ . 

First  determine  the  plug  velocity 
Zone  II 

At<y<  yi 

vlpiug=A^+]-Mryz(0)) 

Zone  IV 

^2  —  y  — 

^W=v'zl,A>%+A2(/t2-^)'2+l 
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Zone  I 


0<y<^ 

vt(y)=KPiug-^-y)Sl+] 

Zone  III 

yt<y<A2 

v}n(y ’)=v"Piug-^2(h-yY2+' 

Zone  V 

A$<  y  <h 


vz  (?)  =  V l! plug  -  A2 (y  -  ^3 )J 


>+l 


Modification  done  for  non-linear  slip  at  the  wall.  The  slip  velocity  at  the  wall  is  given 
by  Us  =  prfz .  Hence,  the  velocities  are  as  follows 

Zone  II 

A,  <  y  <  y(- 

Zone  IV 

A^  <  y<  Ay 

Vz!plug  =  vlplug  +-^2(^2 -3'i)'?2+1 


Zone  I 


0  <  y  <  X\ 

y}iy)=ylpiug-^M-y)Sl+i 

Zone  III 

yi<y<  A 2 

y?l(y)=y"niu8-*2(h-y)S2+' 
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Zone  V 


Aq<  y<h 


v?(y)=vl!Piug-*2(y-A3)S2+' 

Determination  of  volumetric  flow  rate 


On  solving  we  get 


Q 

w 


=  -A. 


Sj  +  2 


htf+] 


- A * 


(^-y/)j2+2+(/»-^r2+2 

S2  +  2 


S-y  +1 


+  U  sfih 


Integration  of  velocity  distribution  for  volumetric  flow  rate 
Volumetric  flow  rate  of  fluid  I  =Qr  +  Qu 

Volumetric  flow  rate  of  fluid  II  =  Q111  +  QIV  +  Qv 
QfiuUL  L  =8/+6;/ 

Q  Fluid  11  5 

i=3 


q!_=vi  ,  MJ|+2 

W  z,plug"\  5)+2 


V/ 


w 


■ = v l  plug  (yi  ) 


Q 


in 


-Vz!plug(M~yi)- 


w 


a  2  (^2 
^2  +  2 


+2 


i/V 


^  -vI!piug(h-‘h) 


w 


g-=yl,Piug{h-M)-^hz^L 

W  z'plug  3/  s2  +  2 


s-j  +2 
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Q  fluid  1  _yll  ^l'V 

W  ~V^yi  ,|+2 


Q 


^r1  =  v"piu8  (h-yt)-  [(h- y  if 2+2 +(h-h) 


5, +2 


Gr 

W 


VL  plug  yi  -  + v (a  -  *• )  -  ^  [(h  -  yi  y>+2 +(h- a,  y> +2 } 


Case  II  (no  zone  2) 

The  details  of  the  derivation  are  not  provided,  but  are  similar  to  case  I,  which  were  given 
previously. 

Condition 


yt  ^  4 


yi 


h-h 

Zone  I 


0  <y<yt 

v!(y)=v"Piug  -a2(^2 -yt)S2+]  -a,U -y)s'+]  -U, -yiY'+]} 

Zone  IV 

Aq  <  y  <  / 13 

vz!plug  ~Pl {Tyz(h)yb2  +  A2(h-  A3Y2+' 

Zone  III 

yt<y<h. 

v!"(y)=vl!plUg-A2(A/2-yy^ 

Zone  V 

VZ  (y) =  Pi (' Tyz  (hjfbl  +  A2  [( h  -  X3 Y2  +l  - (y -  /?3 )52  +l  ] 
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Overall  flow  rate 


^  =  +h(h-j!f,+2jk-W+1 

Case  III  (no  zone  I) 

Condition 

h-h 

Xx  <0 

Determination  of  ryz  (0) 

(^1))S&2  +  A2(/i  —  )‘r2+1  =/?](— r>>z  (0)^' 

Zone  II 

0<y<yt 

Vlplug  =AM0  )^' 

Zone  III 

yi<y</l2 

v{'  ( y) = (-  »*  (0)^ + a2L(4  -  -  W2  -  y)^1  J 

Zone  IV 

Xui<  y  <  X^ 

y  =  X2,vtV  =Vnpiug  =&{Tyz(0jfb'  +A2(X2-yiyi+] 


11 


Zone  V 


A3  <  y  <  h 

v?(y)=v"Piug-*2(y-^Y2+' 


Overall  flow  rate 

+  hft{-Tyz( 0)) 

Case  IV  (no  zone  3) 

Condition 


Qr  _ 

w 


=  A' 


52+2 


yi-^2-  hs 

Determination  of  r,,z(0) 

AM*))*4*  +A2(/r-^)^+l  =A{-Tyz(0)Yb'  +A,(^r'+I 

Zone  I 


0  <  y  <  A] 


V/  GO  =  V/t  p[Ug  -  A  i  -  y  )51  +1  . 

Zone  II 

A,  <  y  <  yt- 

^Lplug  -  A]^' +1  —  P\Tyz  (0) 

Zone  IV 


At  y  =  y,- 

V11  i  =V7  / 

vz,plug  vz,plug 
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At  yi<y<  A 3 

v"  =v"Pi»s=M-*n<f>jfb'  +A,^+I 

Zone  V 

Aq  <  y  <  h 

v?(y)=v;!plug-A2(y-AqY>+' 

Overall  flow  rate 

—  =  U  (- *yz (°)^'  +  +'  ]^ -  Al  aJ^^2*2 

W  ^lv  yz  '  11  1  Si +  2  s2  +  2 

Case  V  (no  zones  1  or  3) 

Condition 

yi<Aq<h 

A,  <0 

Determination  of  7^(0) 

&(»*(*»*■ + a  2(*- V’*'  =[-M^(o)f' 

Zone  II 

OZyZn 

v"=v"Piug  =vlPiug  =[-fryz(0)}b' 

Zone  IV 

yt<y<Aq 

v!V  =  v"plug  =  vlplug  =  [-  PSyz  («)]* 
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Zone  V 


At  y  =  / l3 

Vz  =v"plug 
At  y  =  h 

Vz  =  PliryM^ 

Overall  flow  rate 

^r  =  \ft(-Tyz(0))sb'\i-A2 

Case  VI  (no  zones  1 ,  3,  or  5) 
Condition 

A,  <0 

h-h 

Determination  of  ryz  (0) 

Zone  II 


s2  +  2 


0<y<yi 

vj1  (  y)  =  fi\  (-  Vyz  (fy  fb' 

Zone  IV 

At  y  =  yi 


V. 


II 


z,plug  Vz ,  plug 


I 
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At  yj<y<h 

vlv  =A(-^(°))!t‘ 

Overall  flow  rate 

Case  VII  (no  zone  5) 

Condition 

o<^<yi 

yi<Aq_<h 

A-2>h 

Zone  I 

0  <  _y  <  A, 

V/  =A,U*'+1 

Zone  II 

^  <  3;  <  y(- 

^  =V/,^  =-^,(^(0)^'  +A,^I+1 

Zone  III 

yi<y<Az 

v!"  =A]$+' -M'yzWf'  +a2  [(M-y^'-^-yY1*' 

Zone  IV 

A2<  y</ij 

Vz,v  =  A,^'+1  +A2(^-3'<)i2+1 
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Determination  of  rJZ(0) 

v"plug=fa(±*yz(h)yb2=v"rlug 

Overall  flow  rate 


£-  ■ =  A2[(4  -  y,  Y’- *'  (h  -  y, )- -  ^  "2 1  ■ -  fl  (- r„  (0))* 1, + . A|  L?' *' 1  -  ^ 


J2  +2 


Case  VIII  (no  zone  2  or  5) 

Condition 

4  >  yi 

yt<A^<h 

Aa,>h 

Zone  I 

0  <y<yL 

v[  =^(ryz(0)f'  +A,UI+1  -(A,  -,vf'+l. 

Zone  III 

yi<y<A2 

Zone  IV 

Xi<y<h 

V}v  =02{*yz(°)fbl  -Ai(h-y)S2+] 
Determination  of  t>>z(0) 

Vllplug  =M±hZ(h))Sbl  =Vz,Plug 
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Overall  flow  rate 


<2l  =  a2 

w  z 


(jj zlil 

s2  +  2 


St  +2 


+  /?i(-r 


yz 


(0 )Y>2  3V+A] 


(^-^•ri+2-^i+2 
+  2 


-  A ) 


Case  IX  (no  zones  1  or  5) 

Condition 

4  <0 

yt  <&i 

^<h 

Zone  II 

0<y<yi 

Vlplug  =ft{-*yZ(0)'fbl 

Zone  III 

vl'<.y)=M-*n<.0)'f’'  +  A2L(^ -y,Y‘' -(4  -y)‘»+l  J 

Zone  IV 

y  <  h 

Aty  =  Vf  =  Vl[plug  =  A2(^  -  >^+1  -/?,  (T,z(0)f  • 

Determination  of  7^(0) 

Pi (±  Tyz (h)fb2  =  P\  (-  Tyz (0)Y’'  +  A2  (Aq  -  yi  f1  +l 
Overall  flow  rate 

+fo(±  T^Wf^h-y,) 


£-  =  M-*yzWyb'W-A2 


s2  +  2 


,s-j  +2 
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Case  X  (no  zones  2  or  3) 
Condition 


yt  <  A, 

2 

yi<A3<h 

Zone  I 

0  <y<yiA1<yi 

v!  =M-Tyz(0)Yb'  -A,Lw  -y)I|+l 

Zone  IV 

v/v  =  A2[(h-A3y^  -(y-^Y'+'l+MtyzWf2 

Zone  V 

Ai,  <  y  <h 

Vz  ~ 02(ryz(hjfb2  +  A2L(^-^3)i|+1  -(y~^)'r2+1J 

Determination  of  ryz(0) 

-a,L(^,  -  yi)I,+l -^'■",J=^2(j-jbs(*))’61  +A2L(A-A,)f=H-'. 

Overall  flow  rate 


{h-Ajf^ih-yi)- 


(/t-l3)^+2" 
(^2  +  2) 


+  fii  (-  Tyz  (°  jfbl  yi  +A| 
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Case  XI  (no  zones  2,  3,  or  5) 

Condition 

*  4 x 
'h 

<  3',- 

Zone  1 

0<y<yi;/L2<yi 

vl  = /?,(-*„(<>  )fi  -a,Lu,  -,)•<+' 

Zone  IV 

yi<y<X3 

At. }  >  h 
h  <y  -i 

v"  = /hKm'f' 

Determination  of  tyz(0) 

M-^n(o)'fb<  -a,L« =Mt„(h)Y^ 

Overall  flow  rate 


-  P\  (~  Tyz  (0)f>  yi  +  A] 


~(4-y,-)J|+2-ff+2~ 

5]  4*  2 


Case  XII  (no  zones  3  or  5) 
Condition 

/^>h 
h.  <  yi 
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Zone  1 


0  <y<A\ 

v!  =-a,Lm  -?)'■*' 

Zone  II 

A\<y<  yt 

V"  =vlplug  =A(-^(0)f'  +A1^I+I 

Zone  IV 

y" = y !! Piug  =  M* +1  -  p\  {ryz  (o)f' 

Determination  of  t>iz(0) 

y l plug  =  P\  (“  Zyz  (°))r/;'  +Al^'+I  =/^2 {Tyz(h)ybl  =Vz,plug 

Overall  flow  rate 

2L=h/i2(±ty,wYb>  -  a,  -3— 

W  1  5, +2 

Conditions  to  determine  which  case  applies 
Case  I  (all  zones) 

0<A,  <yh  yi<^i,  Ay  <h 
Case  II  (no  zone  2) 

yt<  A]t  y;<X 2,  Ay  <h 
Case  III  (no  zone  1) 

yi  <  Aq  ,  Ay  <  h ,  Ay  <  0 

Case  IV  (no  zone  3) 

0</ti<^,  yi<A2<h,  >?2  <  yi 
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Case  V  (no  zone  1  or  3) 

yi  <  Ay  <  h ,  Aq  <  y;,  X\  <  0 
Case  VI  (no  zone  1 ,  3,  or  5) 

^  0 ,  —  yi  >  Ay  —  h 

Case  VII  (no  zone  5) 

0  <  ^  <  yi ,  yi<  Xi<h  ,  Ay>h 
Case  VIII  (no  zone  2  or  5) 

A\  > y,- .  yi^A^<h,  >h 
Case  IX  (no  zone  1  or  5) 

^<0,  yi  <  A2,  Ay  <  h 
Case  X  (no  zones  2  or  3) 

yt<A\,  yi>A2,  yt  <  Ay  <h 
Case  XI  (no  zones  2,  3,  or  5) 
y^Aj,  Ay>h,  X1<  yi 
Case  XII  (no  zones  3  or  5) 

0<A\  <yit  Ay  >  h  ,  ArL<yi 


Typical  comparisons  of  the  experimental  results  versus  the  predictions  of  the  model  are 
shown  in  figures  31  to  34.  Overall,  the  predictions  of  the  flow  rate  versus  the  pressure  drop  are 
adequate. 


SUMMARY  OF  GENERAL  CONCLUSIONS  AND  LESSONS 

•  The  manufacturing  of  co-extruded  grains  of  highly  filled  propellants  is  an  intricate 
process,  which  requires  a  realistic  understanding  of  the  various  types  of  interface, 
surface  and  bulk  instabilities. 

•  The  interface  between  the  two  fluids  can  shift  during  the  co-extrusion  process  if  the 
wall  slip  and  the  shear  viscosity  of  the  two  phases  are  not  similar. 

•  One  or  both  of  the  two  surfaces  of  the  co-extruded  grains  can  be  subjected  to  flow 
instabilities  and  the  related  development  of  gross  surface  irregularities,  which  can 
render  the  products  unacceptable. 
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•  The  use  of  mathematical  modeling  allows  the  understanding  of  the  details  of  the 
process  to  make  educated  designs  and  selection  of  processing  conditions  possible, 
however,  these  models  need  to  be  further  expanded,  validated  and  converted  into 
tools  that  are  user-friendly  for  the  day  to  day  workings  of  the  Department  of 
Defense. 

Specific  Conclusions  and  Lessons 

Material  Characterization 

•  Consistent  with  what  we  have  advocated  since  1986,  the  wall-slip  behavior 
and  the  viscoplasticity  of  the  phases  to  be  co-extruded  need  to  be  accurately 
determined.  It  is  now  clear  that  the  platform  that  we  have  established  is 
equally  valid  for  the  co-extrusion  process. 

•  The  critical  shear  stress  at  which  the  transition  from  stick  to  slip  for  a  given 
suspension  should  be  an  additional  parameter  to  be  characterized  during  the 
characterization  of  the  phases. 

•  The  compressibility  of  the  melt/suspension  and  the  pressure  dependence  of 
the  slip  coefficients  need  to  be  incorporated  into  the  mathematical  models  of 
the  co-extrusion  process  to  accurately  determine  the  conditions  under  which 
flow  instabilities  are  on-set,  so  that  such  instabilities  can  be  minimized  or 
eliminated. 

Formulation  for  Processability 

•  The  formulations  of  the  two  phases  to  be  co-extruded  should  be  tailored  if 
possible  to  match  the  wall  slip  behavior  and  the  shear  viscosities  of  the  two 
phases  to  be  co-extruded. 

•  If  the  differences  in  the  rheological  behavior  of  the  two  phases  are  not 
significant  then  the  interface  remains  stable  during  the  co-extrusion  flow. 

•  If  the  wall-slip  condition  remains  stable  during  processing  flow  instabilities  are 
not  encountered  and  the  formulation  affects  the  stability  of  the  wall  slip 
condition  also. 

Mathematical  Modeling 

•  We  have  derived  a  complete  set  of  equations  to  describe  the  co-extrusion  of 
viscoplastic  fluids  subject  to  wall  slip  in  various  configurations,  including  the 
side  by  side  co-extrusion  of  two  viscoplastic  fluids,  a  case  studied  extensively 
with  our  experimental  setup. 

•  The  comparisons  of  the  results  of  the  experimental  studies  with  the 
mathematical  modeling  results  is  acceptable  at  moderate  rates,  however, 
problems  are  noted  at  the  low  and  high  flow  rate  regions. 
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•  The  sources  of  the  discrepancies  for  co-extrusion  at  the  low  and  high  flow 
rates  need  to  be  determined,  especially  by  revisiting  the  fits  of  the  wall  slip  and 
shear  viscosity  characterizations  at  the  low  and  high  shear  rates  (parameters 
generally  do  not  cover  the  entire  range  of  shear  rates). 

•  We  have  developed  a  new  model  to  predict  the  on-set  of  flow  instabilities  on 
the  basis  of  the  compressibility  of  the  melt/suspension  and  the  change  in  the 
wall  slip  coefficient  with  pressure  along  the  length  of  the  die,  which  appears  to 
be  very  accurate  for  extrusion  flows. 

•  However,  the  compressible  fluid  model  needs  to  be  expanded  into  the  co¬ 
extrusion  process. 

•  At  the  first  opportunity  possible,  the  steady  FEM  models  of  SIT  will  need  to  be 
converted  into  time  dependent  models  to  allow  time  dependent  calculations  to 
be  made  for  process  control  purposes  (model  based  process  control)  and 
better  design  of  the  process.  Such  calculations  can  also  be  used  to  generate 
functionally  graded  co-extruded  propellants  in  the  axial  direction  in  a  cyclic 
fashion. 

Design  and  Operating  Condition  Selections 

•  If  the  formulations  cannot  be  tailored  to  match  the  rheological  behavior  then 
the  geometries  and  operating  conditions  should  be  tailored  to  allow  the 
ultimate  rheological  behavior  of  the  two  phases  in  the  co-extrusion  die  to 
match.  This  can  be  accomplished  to  some  extent,  for  example,  by  the 
selection  of  two  different  temperatures  for  the  two  phases. 

•  It  is  the  presence  of  the  converging  section  in  the  design  that  we  have 
selected  which  has  given  rise  to  the  changing  interface  location  as  a  function 
of  time  during  co-extrusion.  More  of  PDMS  has  flown  in  the  die  and  thus  more 
of  it  is  depleted  from  the  reservoir  section  during  the  course  of  the  co-extrusion 
followed  by  the  depletion  of  the  pure  PDMS  and  then  the  consequent  increase 
in  the  flow  rate  of  the  filled  PDMS  which  now  occupies  more  of  the  converging 
section  of  the  die. 

•  A  reservoir  region  adjacent  to  the  fully  developed  section  of  the  die  should 
never  be  included  in  the  design  of  the  co-extrusion  die. 

•  One  needs  to  determine  apriori  flow  conditions,  which  provide  extrudates  that 
are  free  of  surface  irregularities  on  their  entire  perimeter. 

On-line  Process  and  Product  Quality  Control 

•  The  pressure  readings  change  very  little  during  the  process,  however,  the 
change  in  the  pressure  gradient  is  significant  along  with  a  significant  change  in 
the  interface  location. 
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•  Thus,  during  the  co-extrusion  process  one  either  needs  to  determine  the 
pressure  values  very  accurately  to  generate  accurate  pressure  gradient  values 
and/or  continuously  follow  the  interface  location  as  a  function  of  time  during 
the  co-extrusion  process. 

•  Currently,  there  is  no  method  to  follow  on  an  on-line  basis  the  interface 
location  during  co-extrusion.  This  problem  needs  to  be  addressed  as  a  quality 
control  issue. 

'  •  If  the  temperature  distributions  of  the  two  streams  change  as  a  function  of  time 
during  the  co-extrusion  process  the  interface  location  is  likely  to  change  also, 
since  the  shear  viscosity  and  wall  slip  will  also  change.  Thus,  the  operating 
conditions  need  to  be  carefully  selected  and  then  maintained  during  the 
process. 
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Figure  1 

Principal  experimental  arrangement  for  approaching  the  co-extrusion  problem;  off-line  slit 

rheometer 


Cartridge  with  the 
two  suspensions  placed 
side  by  side  as  if  they 
are  fed  by  two  extruders 


Pressure 

transducers 


J 


Figure  2 

Apparatus  for  the  investigation  of  the  shape  of  the  interface  and  flow/interface  instabilities 
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Figure  3 

Placement  of  a  partition  in  the  mid-plane  of  the  cartridge  and  the  two  sides  for  allowing  the 

interface  to  be  traced 


Figure  4 

Partition  and  cartridge  used  in  the  experiments 


26 


Suspension  of  10-40°/o  glass 

Pure  PDMS 

In  white  color 

Translucent 

\ 

/ 

< 

\v 

1] 

Figure  5 

Coloring  of  the  two  different  co-extrusion  layers 


Figure  6 

Exit  of  the  co-extruded  grain  from  the  slit  die 


Results 

Pure  PDMS  and  10%  Hollow  Glass 
Filled  PDMS 

Co-Extrusion 


Figure  7 

Results  for  0  and  10%  co-extrusion 


Figure  8 

Pressure  versus  distance  in  the  slit  die  for  different  flow  rates  (cross-head  speeds) 
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t; 


Pure  PDMS  and  10%Glass+90%PDMS  Co-Extrusion  in  the  offline  Slit  die 

Gap=7.05,  mm  T~23°C 


Figure  9 

Distribution  of  the  two-co-extruded  layers  during  extrusion:  0.1  in./min,  cut  in  the  transverse  to 

flow  direction 


No  extrudate  distortions:  Wall  shear  stress  around  6700  Pa. 


10%  Hollow  Glass+90%  PDMS 


Pure  PDMS  and  10%Glass+90%PDMS  Co-Extrusion  in  the  offline  Slit  die 


Figure  1 0 

Distribution  of  the  two  co-extruded  layers  during  extrusion:  0.1  in./min,  cut  in  the  transverse  to 

flow  direction 
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Figure  1 1 

Distribution  of  the  two  co-extruded  layers  during  extrusion:  0.1  in. /min,  cut  in  the  parallel  to  flow 

direction 


For  pure  PDMS  and  10% 
Glass  Filled  PDMS 


•No  change  in  the  interface  location  with  time  during  coextrusion 

•The  flow  instabilities  are  onset  around  60,000  Pa, which  agree 
with  the  behavior  of  the  individual  materials. 


Figure  12 

Summary  of  co-extrusion  results  for  pure  PDMS  and  10%  by  volume  glass 
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Results 

Pure  PDMS  and  40%  by  volume 
spherical  glass  filled  PDMS 
Co-Extrusion 


Figure  13 

Results  for  0  and  40% 
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Figure  14 

Pressure  versus  time  for  co-extrusion  of  0  and  40% 
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Temperature,  °C 
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Figure  15 

Pressure  versus  distance  in  the  slit  die  for  different  flow  rates  (cross-head  speeds)  during  co¬ 
extrusion  of  0  and  40% 


Figure  16 

Distribution  of  the  two  co-extruded  layers  during  extrusion:  0.1  in. /min,  cut  in  the  transverse  to 

flow  direction 
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Figure  17 

Distribution  of  the  two  co-extruded  layers  during  extrusion:  0.1  in. /min,  cut  in  the  parallel  to  flow 

direction 


For  pure  PDMS  and  40% 
Glass  Filled  PDMS 

•A  shift  in  the  interface  location  with  time  is  observed  for  all 
conditions. 

•The  flow  instabilities  are  onset  around  60,000  Pa  for  pure 
PDMS  and  only  trace  amount  of  surface  irregularities  are 
observed  for  40%  filled  PDMS  agreeing  with  the 
behavior  of  the  pure  individual  components 


Figure  18 

Summary  of  co-extrusion  results  for  co-extrusion  of  0  and  40% 
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y=  h,  gap  of  the  slit  die 


Fluid  1 


y=0 

Schematic  of  co-extrusion  of  two  fluids  flowing  side  by  side. 

Figure  19 

For  the  mathematical  modeling:  basic  flow  configuration  for  co-extrusion 

PDMS  with  40%  glass  spheres: 

Slip  velocity.  Us  vs  wall  shear  stress,  xw 

us  =PTwsb(1  +  tanh(a(xw  -Ta))) 


Shear  stress.  Pa 


Figure  20 

Slip  velocity  versus  wall  shear  stress  for  PDMS  binder  with  40%  by  volume  filler 
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Two  viscoplastic  fluids  subject  to  wall  slip-  Case-I 


Figure  21 

Shear  stress  distribution  for  case  I:  co-extrusion  of  two  viscoplastic  suspensions  subject  to  wall 

slip 
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Figure  22 

Velocity  and  shear  stress  distributions  for  case  I:  co-extrusion  of  two  viscoplastic  suspensions 

subject  to  wall  slip 
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Figure  23 

Velocity  and  shear  stress  distributions  for  case  IV:  co-extrusion  of  two  viscoplastic  suspensions 

subject  to  wall  slip 
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Figure  24 

Velocity  and  shear  stress  distributions  for  case  VII:  co-extrusion  of  two  viscoplastic  suspensions 

subject  to  wall  slip 
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Decision  tree  for  case  determination  for  the  analytical  solution  of  the  co-extrusion  flow 


Predictions  aiming  at  the 
experimental  conditions  for 
the  co-extrusion  of  pure  PDMS 
and  40%  glass 
a.  Interface  at  the  mid-point 

Figure  26 

Predictions  for  co-extrusion  of  pure  PDMS  binder  with  PDMS  with  40%  filler 
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Figure  27 

Predictions  of  velocity  and  shear  stress  distributions  for  co-extrusion  at  ram  velocity  of  0.1 

in. /min 


Cross-head  speed=l  in/min.  dP/dz.^ — -12  MPa/m,  QprtdJctio„=  5.14*10  7  m3/s 

Figure  28 

Predictions  of  velocity  and  shear  stress  distributions  for  co-extrusion  at  ram  velocity  of  1  in. /min 
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Figure  29 

Predictions  of  velocity  and  shear  stress  distributions  for  co-extrusion  at  ram  velocity  of  5  in. /min 


Summary  of  experimental  data  and  predictions 
for  pure  PDMS  and  40%  Glass  Filled  PDMS 


Crosshead  speed 

Pressure 

Flow  rate 

Flowrate 

of  the  ram  extruder 

Gradient, 

actual. 

predicted. 

inches /min 

MPa/m 

m3/s 

m3/s 

0.1 

-2.9 

5.40E-08 

1.80E-08 

1 

-12 

5.40E-07 

5.14E-07 

5 

-19 

2.71E-06 

2.71E-06 

10(a) 

-23 

5.4GE-06 

4.00E-06 

10(b) 

-26.6 

5.4GE-06 

5.42E-06 

Figure  30 

Experiments  versus  co-extrusion  theory 
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Predictions  aiming  at  the 
experimental  conditions  for 
the  co-extrusion  of  pure  PDMS 
and  40%  glass 

b.  Interface  between  the  two 
fluids  changing  with  time 

Figure  31 
Interface  location 


Figure  32 

Pressure  gradient  versus  time  at  5  in. /min 
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Figure  33 

Velocity  and  shear  stress  profiles  at  5  in. /min 


Figure  34 

Experimental  versus  theory:  pressure  versus  time  in  the  slip  die  during  co-extrusion 
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